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Abstract 

The theory for manifestation of stochastic appearance of embryos 
in the global decay of metastable phase has been constructed. The 
regime of droplets growth is supposed to be both free molecular one 
and diffusion one. The deviation for a mean droplets number from the 
value predicted by the theory based on averaged characteristics have 
been calculated. The value of dispersion for the distribution of the 
total droplets number in particular attempt has been also calculated 
analytically. Comparison with results of numerical simulation has 
been given and the correspondence between simulation and analytical 
approximate results is rather good. 

The aim of this paper is to give analysis of decay in diffusion regime of 
the droplets growth in frames of the method given in [!Q and corrected in [H] . 
This method is based on the following approximation 

• Until the half of the nucleation period (in the time scale) the supersat- 
uration remains unperturbed and the droplets appear independently. 

• In the second half of the nucleation period the droplets appear under 
the supersaturation changed by the vapor consumption by the droplets 
appeared in the first part of the nucleation period. 

The part of the droplets size spectrum which corresponds to the droplets 
appeared during the first part of the nucleation period will be called as the 
first part of spectrum, the remaining part will be the second part of spectrum. 
Here this model will be generalized by means of the following changes: 
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• The duration of the first part of nucleation period will be not exactly 
the first half but it will be regulated by some parameter p instead of 
1/2. 

• The regime of growth will be an arbitrary power law of growth (the 
power type is adopted only to present concrete formulas, it isn't prin- 
cipal here). 

• It is supposed that the distribution of droplets in the first part of spec- 
trum can differ from the distribution 

P (iV)~exp( l 2 - > ) 

of N independently appeared droplets during the period of time corre- 
sponding to the mean value N of appeared droplets. 

Now instead of the previous distribution the distribution with renor- 
malized dispersion 

will be used. Here i[> is parameter of renomalization. This change is 
initiated by the using of the property of self similarity of spectrums in 
the first iteration (see [S]) used for investigation of stochastic effects in 

The structure of consideration will be the following 

• At first we shall present the formal generalized model 

• Then the parameter p of the boundary between the first and the second 
parts will be calculated. 

• In the last part the effect of the "growing volumes of interaction" will 
be described and it will be shown how to include this effect into already 
presented scheme. 

• In Appendix the theory for the free molecular regime of growth will 
be presented. This theory was given in [Hj but now some modifications 
and improvements have been made. So, a new variant can be found in 
appendix. 
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The main object of investigation will be the calculation of dispersion 
of the droplets distribution since the mean value will be very close to the 
characteristic calculated in frames of the theory based on the averaged char- 
acteristics (TAC), all arguments given for the free molecular regime of the 
droplets growth remain here practically the same the only principally 
new feature is to use the theory with explicit profiles of density around the 
growing droplets 

The smallness of fluctuations which will be shown here plays a very im- 
portant role in transition from the free molecular kinetics to the diffusion 
kinetics. Really, such transition was shown in [3] only in general features on 
the level of TAC and the possible giant fluctuations initiated by stochastic 
appearance of droplets would destroy this transition in frames of the known 
approaches. Fortunately the giant fluctuations don't appear and the justi- 
fication of their absence lies in the moderate value of dispersion calculated 
here. 

1 Calculation of dispersion 

Here we shall analyze the situation with diffusion regime of metastable phase 
consumption. At it has been already noticed this situation is characterized 
by the growing volumes of interaction which produces additional difficulties. 
The growing volume of interaction means that the fixed point in the volume 
will be under influence of vapor consumers which appeared in a volume V 
with a center at this point and with a radius y/ADt where t is a time counted 
from the consumer appearance up to a current moment and D is a diffu- 
sion coefficient. Certainly the last value can be regarded only as the rough 
estimate, the precise expression can be obtained only of the base of precise 
profiles of a substance gap [I]. 

At first we briefly follow the way to derive the estimate for dispersion 
proposed in l a to use it for diffusion regime of metastable phase consumption. 
As we have noticed the formalism of a volume of interaction is a rather rough 
estimate (precise theory has to be based on a form of substance gap from 0j). 
Thus, the accuracy here is limited. So, in frames of adopted accuracy we shall 
use here only some rough way to estimate dispersion of a distribution and 
can use formalism from £Q with appropriate comments and modifications. 
This way is useful because it allows to understand formulas from £Q which 
were widely used in above consideration. 
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Here we suppose that the rate of embryos growth is given by 



dv 
~dt 



where v is the number of molecules inside the embryo and s is a power 
parameter. We suppose that s takes values which are essentially greater 
than (s isn't too small in comparison with 1) and essentially less than 1 
(1 — s isn't too small in comparison with 1). For the free molecular regime 
of metastable phase consumption s = 2/3, for the diffusion regime s = 1/3. 
Then the characteristic p of a droplet defined as 

p = v^ +1 

grows with velocity independent on it's value 

dp/dt ~ 1 

The period of nucleation is divided in many elementary intervals with a 
length A (or Aj) in a p-scale, the average total number of interval in the 
whole nucleation period is M, the number of embryos iVj appeared during 
every elementary interval is supposed to be big iV, 3> 1. 

The fact that the total number of intervals is M and their length is A (in 
estimates we can put it one and the same for all intervals, at least we can 
take the smallest length and attribute it to all intervals) means that 

r/(MA)^ +1 > ~ (-J- + i)c(o) , 

1 — s 

where ((0) is the initial supersaturation, / is the initial average amplitude 
of embryos size spectrum [3], F is parameter from TAC Namely the last 
relation can be derived from the first iteration in TAC. It will be used to 
express / through M, A. 

During every interval the increase Sp of the value p of already existing 
droplet will be 

5p ~ AC(0)/r 

with some characteristic time r which will be dropped out in final formulas. 

Denote by iVj the number of embryos appeared at i-th interval. The value 
Ni is the average iVj, Pi(Ni) is the density distribution on N { . 
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Denote by N^ % ' the number of embryos appeared at the first i-th intervals 
A. The value N^' is the average N^ l \ pW(A^W) is the density distribution 
of N®. It is important to stress that isn't equal to N® which is the 
corresponding value completely calculated in frames of TAC (stochastically 
appeared droplets in the previous moments of time taken into account). 

At the first P intervals (i.e. at the first part of nucleation period) the 
metastable substance consumption is negligible. Here P = pM and p is 
parameter of constructions. Then for k < P the distribution P^ of the 
number of droplets N^ k > appeared in the first k intervals is 

«~-(-^^> 

and it is is a standard Gaussian distribution of k independent subsystems 1 . 

We introduce here the parameter of renormalization if> to be able to fulfill 
calculations suggested in [Hj based on the approximate similarity of nucleation 
conditions 

For k > P the expression for P^ k > will be formally the following 

pW(iv«) = £ {[P^N^n 

Ni, ,N P i=l 

where 

P(iV J )~exp(- 



2^Ni ' 

k — 1 k—1 

n= e n pjmp^-j:^) 

N P+1 ,....,N k _ l j=P+l 1=1 

Here the function fl is extracted because it corresponds to the second part of 
spectrum. The function fl is superposition of independent subsystems and 
can be calculated explicitly 



(N^-ZliNt-EUp+iNj] 
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O / \ '1=1 ' *—'j = ~~t1 J' \ 

\l ~ exp r =?■ 

* 2^j=p+i iv j 

The problem is that the values Nj standing in the last expression are 
unknown now. They are functions of all Ni, i = 1,..,P (since the two 

^^Here and later without any notation we don't keep the preexponential factors before 
Gaussian distributions, they can be easily reconstructed by integration. 
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parts model is used). So, the distributions from the first group i < P have 
the influence on Q. This influence is given through the values Nj. Now 
we shall get an expression for this value. At first we recall the exponential 
approximation which can be found elsewhere (see [3]) 

so 

Here Gj is the number of molecules in a liquid phase at j-th interval. For 
this value we can write 

i=l 

where pf' is the " size" p of a droplet appeared at z-th interval at the moment 
corresponding to the j-th interval. Then 

Gj = X>i((ACo/T)j " (ACo/r)*)^ = X>.(j " i)^(ACo/r)^ 
i=i i=i 

and 

Nj = ^ 1 exp(-^(ACo/r)^EiVi(i-i) T ^) 

i= i 

Now we have to get another expression for 

K=^(ACo/t)t^ 
so 

which appeared in the last relation. One has to clarify the meaning of M 
as the total number of intervals. Certainly, it is no more than a conditional 
definition to say that the end of nucleation is the moment when the averaged 
rate of nucleation falls e times in comparison with initial value. Then the 
corresponding equation will be 

—Gm = 1 
so 

We have to calculate Gm in some approximation. Let it be the ideal approx- 
imation, i.e. the first iteration [3.,. Then 

M 

KY,Ni(M -i) l/{l - s) = 1 

i=l 
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or 

M 



KN^M 

For M ^> 1 summation can be replaced by integration which gives 

M /-M ]Vf 1/(1 " s)+1 



and 



Then 



Here 



Y,{M - i) 1/{1 - s) = {M -xf'^dx 
i=i ^° ~ + 



*=(*i _l r i r 1 (1) 
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i=l 



iViM ( ^ +1) 

is simply the scaling factor. One can argue whether the style of calculation 
of Gm is appropriate. Really we used another approximation which differs 
from the two cycle model which we are going to use. But there is absolutely 
no sense how to calculate the last value. One can simply say that our model 
is to take the first part as P = pM where the definition of M is given by (JTJ. 
The problem is how to choose p. 

We continue to simplify expression for Nj, j > P. We recall that the 
droplets from the first part only consume vapor. Then we come to 

Nj = N x exp(- * £ Ni-*— + - i)^=) (2) 
iV\M ( ~ +1) t=x l ~ s 

Now we shall simplify the last formula. We make the following transfor- 
mations: 

- we add and subtract Ni, which is the ideal mean value of droplets 
appeared during the elementary interval. Then 

ft - ft «p(-^£ ( i + %^)(^- s + m - .•)*) 
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and 

p 



We can calculate £*=i + l)(j - i)^. For M > 1 

p 1 1 r p 

E(i + l)(j = (i + 1 )/ (j'-a;)^tfa; 

~i L — s I — s Jo 

and since we suppose 1/3 < s < 2/3 we have 1/(1 — s) ^> 1 and 

/ (j - x)^dx = (j - x)~^dx = (- h l) _1 j i^ +1 

jo jo 1 — s 

for j which aren't too big in comparison with P (namely the parameter p 
will be chosen in such a way that P can not be too small in comparison with 
M and certainly j can not be essentially greater than M). 
Then 

1 I i 

and 

Now we linearize the last exponent over 

This linearization is possible and it differs from linearization over 

1 Ni-N lf 1 



M ( ^ +1) iVi 4-s 



+ i)0'-») T = 



made in PQ. The linearization in PQ since we study fluctuations namely in iVj 
is very doubtful. In our approach we can see the compensation of fluctuation 
during the period comparable with the whole nucleation period and, thus, 
the linearization is possible. 

It is necessary to stress that after the linearization one can not pretend 
on the derivation of deviation of the mean value of droplets from the value 
prescribed by TAC. Contrary to P we don't pretend on this but here we 
are interested only in dispersion. To prove the smallness of deviation of the 
mean value of total droplets from the value given by TAC one can use other 
approaches (see 0. 

Here we directly use (j^- + l)(j — i)~ as the quantity of substance in 
droplets of a given size (the relative deviation of the size during the ele- 
mentary interval is small) instead of inappropriate integration over constant 
distribution (it is stochastic and isn't constant) in pQ. 

In the zero approximation we get 

In the first approximation 

1 Ni — N-i 1 i 

Nj = N, eM-i-r—, )(1 - — — £ "V^ 1 + 1 W ~ 0^ 

M — + 1 M~ +1 ~i ^ 1-s 

Now we know the influence of embryos of big sizes and can directly analyze 
the numerator of Gaussian distribution for Q. We have 

P k P P k 

(AT«-£AT.- £ N,) = {N^-Y J {N l -N l )-Y,N l - E Nj) 

i=l j=P+l i=l i=l j=P+l 

P P k 

t=i i=i j=p+i 

P k 

^ N (k)_^ N ._ Ni)+pNi _ Nj) 

i=l j=P+l 

Now we shall derive formula for Yjj=p+\ Nj. We start with (J2J) where N = iVi 
and get 

E ^ = ^o E ^(-^iE^ + iifi-.) 1 ^) 

j=p+i j=p+i MM ; iV o i=i l-s 
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Since M > 1, P ~ M 
p 



jr(— + i)(j - iy/(i-)/ M w-s)+i ~ jW-8)+i/ M i/(i-s)+i 
i=i i — s 



So, one can come to 



E *i = *o E exp(-( J 7M)- 1/(1 - s)+ >xp(- ¥T7 i^ 
j=p+i j=p+i JW 



1 £(-!-- + l)(j - ,) 1/(1 " s) iV l ) exp(£ (J- + l)(j - i)W-.)l M W-^) 
iV o i=1 1 - s i=1 1 - s 



Then 

E = E exp(-( J /M) Ml/(1 - s)+1 ) 
j=p+i j=P+i 

2=1 u 



Now one can fulfill linearization of exponent over 

p 

which leads to 

k k 



i yv J_ + i )(7 - _ i/(i-) W-*o) 



E *; = *o E exp(-( J /M) Ml/(1 - s)+1 ) 
j=p+i j=P+i 

The formula is ready. 
Then 

P k P 

(iVW - ^(7V. - 7V0 - PTVi - £ N j ) = (N^-'£(N i -N 1 )-PN 1 - 

i=l i=P+l »=1 

J»„ £ exp(-07M)W-.)«) ( l-^-^f(^ + l)0-i)V( 1 -.)(«) 

j=P+l JW i=l 1 6 iV o 

10 



)) 



The next transformation gives 

(N^-^-Nj-PN.-N, £ exp(-07M) 1 /a-^i)(l- J_ s)+i 

fl^ + ilO-.) 1 '' 1 -''^^))- 

p fc 
- £(jv; - Nl ) - PN, - N £ exp(-07M) 1 /d-)+i) + 
i=i j=p+i 

J»o _E +i exp(-07M)W-.m )]gI7 ^f(-l- + 1)0 - ;,!/(.-., W^M, 
Having noticed that 

AT( fc ) = iV 1 (P+ £ exp(- 7 , — ) 
i=P+ i M ( — +1) 

is the mean number of droplets prescribed in the zero approximation 
by TAC, we come to 

P k 

{N (k)_Y / (N i -N 1 )-PN 1 -N ]T exp(-07M) 1 /a-)+i) + 
i=i i=P+i 
fc 1 p 



j=P+l JW i=l 1 6 iV o 

0V (fe) - - #i) - JVW + 7V E exp(-(i/M) 1 /(i-)^) 

i=i i=P+i 

Now we shall change the order of summation and come to 

p 

i=\ 

N t exp(-07M) 1/ ' 1 -» +I ) S T 7 ^ TT f:( r i 7 + 1)0 " w -> = 

j=P+l JW i=l 1 6 Jv o 
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(iv (fc) -J2( N i~ N i)- N (k) + 

i=l 

- n ) E exp(-07M)V(i-^) Mi/( |_ <)+i (r L- + - ova-); 

The last transformation leads to 

p 

( N (k) _J2(N i -N 1 )- + 

i=l 

- iV ) E eM-(j/Mf^ )+1 ) I _. )+1 (t^- + l)(j " V(1 - S) ) = 
i=i j=p+i 

p 

(JV« - iV (fc) - ^(iV, - jV )(l - 



i=i 



E ex P(-(^ M ) 1/( " s)+1 )]^FT(T^ + W - ^) 1/(1 - s) )) 

J=p+1 



or 



I , 1 



£(iV,-iV ) E expl-OYM)^-^) ( + 1) (J 

i=i j=p+i 6 



-i)V(i-)) 



i=i 

where 



af» = (1 - ± exp(-07A/)W--»«)- gI7 i- y - T (_L- + ,)0 - W -») 



j=P+l 

The last definition can be rewritten as 

«?> = (1 " &?>) 

where 



*f° = . E ex p(-0V M ) 1/(1 " s)+1 )]^v(^(T^ + 1)0" - 



i=p+l 
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Now we shall analyze denominator of Gaussian distribution for f2. Here 
we can use Nj in the zero approximation and come to 

k 

Nj = N (k) - PNx 

j=p+i 

The use of zero approximation is based on the following simple estimate. Re- 
ally, the last linearization is based on the smallness mainly of the parameter 
which can be simply neglected in denominator since 

E exp(-(A)i/(i-.) + i) 
j=p+i m 

isn't too small. 

Now we can directly go to the calculation of pW(N^) for k > P. The 
starting formula is 



p 

P (fe)(7V«) = f°° d{N x - Ni)... f°° d{N P - N P ) f[exp 

J —OO J — OO • i 



exp( 



(N-N,) 2 , 



(3) 



2(JVW - PN-l) 



where we because of N ^> 1 have replaced summation by integration. The 
limits of integration can be infinite ones because of N ^> 1 and estimates 
predicted by Gaussian distribution. 

To fulfill integration we have to use a simple formula 



/oo x % ([x -\- d) 2 d? 
exp(-— ) exp( )dx ~ exp(-— — -) 
-oo C + l z C 

We have to use this formula P times with 

x ~ Xi = N — Ni 



c ~ a = 2ipN l 
L ~ Li = —a] 
b~bi = 2(iV (fe) - PN X ) 
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This leads to 

where dispersion is given by 

£>(*) =iV( fc )-(P-^£af )2 )iV 1 



i=l 



Now we can calculate dispersion at infinite k. Since M > 1 we can 
substitute summation by integration. The number of interval corresponds to 
the variable z * M, then z — 1 corresponds to the end of nucleation. 

Since 

jV(oo) = at 1 (p + £ exp (_(A)i/(i-s)+i)) ~ Mi y ia 
j=p+i M 

a = exp(-x 1/(1 - s)+1 )dx 
Jo 

then 

One can easily calculate 6^°°'' : 

-> (1/(1 - s) + 1) / exp(-a: 1 /( 1 - s ) +1 )(x - yf'^dx 
Also one can note that 

p 



MU Jo 
Then 

£,(oo) = - (1 - - 

where 
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r-p poo 

(h = 2(1/(1 - s) + 1) / dy dxeM-x 1/(1 ~ s)+1 )(x ~ y) 1,(1 ~ s) 

Jo Jp 

f-'P POO 

f3 2 = (1/(1 - s ) + If / dy( / dxexp{-x 1/(1 - s)+1 ){x - y) 1/(1 " s) ) 2 

Jo Jp 

This result can be directly calculated. 

2 Calculations of parameters of the model 

To compare analytical constructions with the real value of dispersion we shall 
present numerical simulations. 

The results for the mean value of droplets calculated in units of the 
droplets number calculated in TAC as function of the droplets number cal- 
culated in TAC are presented in figure 1. 
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Figure 1: Mean value of the droplets number 
The relative dispersion as a function of the mean number of droplets is 



shown in figure 2. 
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Figure 2: Dispersion 

Now we shall perform calculations. We study equation 

5 f z 

9(x) = o / ex P(-9(x))(z - xf' 2 dx 
I Jo 

Namely the coefficient | corresponds to the final formulas in the previous 
section. 

At first we calculate the form of spectrum. It is shown in figure 3. Two 
curves are drawn. The upper one in the second iteration approximation, the 
lower curve is the first iteration approximation. The spectrum lies between 
these curves. 

One can see that the characteristic length of spectrum is x ~ 1.2 4- 1.25. 
The total number of droplets is N TAC = 0.9292. The first iteration gives 
N first = 0.8773. 

Now we investigate the two cycle model. The relative number of droplets 
as a function of the boundary p is shown in figure 4. 
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Figure 3: The form of spectrum in TAC 



Here two curves are drawn. The upper corresponds to the number of 
droplets referred to the number of droplets in the first iteration, the lower 
curve is referred to the precise value of the droplets number. One can see 
that the deviation of the second curve is two times closer than the deviation 
of the first curve in the value of minimum. The minimum is attained at 
Po = 0.64 which is approximately the 55 percent of Xq (the same part as in 
the free molecular regime! |6]). 

Now we calculate (3. The values (3% and fa are drawn in figure 5. One 
can see that the behavior of /3 is mainly governed by behavior of fa. 

The ratio fa/ fa is drawn in figure 6. 

One can see that there is no extremum of the ratio. So, there is no special 
boundary which provides the quickest convergence of the chain fa, fa, ■■■ 
and one has to use po = 0.64 going from the most precise approximation for 
number of droplets as a criterium for such boundary. Also one can see that 
the ratio is small for all p. 
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Figure 4: Mean value of the relative droplets number in two cycle model 

The value of relative dispersion calculated according to 

D = 1-13 /a 

is drawn in figure 7. 

There are two curves: the upper is dispersion calculated with precise 
number of droplets, the lower is dispersion calculated with the number of 
droplets calculated in the first iteration. Both curves have extremum at 
p = 0.68. Even for the lowest curve all values are greater than 0.54 while the 
real value is 0.44. 

If we take into account that a also has to be calculated in two cycle 
model. Then the value of relative dispersion 2 as a function of the boundary 
p ini the two cycle model is drawn in figure 8. 

One can see that practically nothing has been changed. Minimum is at 
0.67, the value of minimum is near 0.6. 

2 More precisely there is a square of relative dispersion. 
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Figure 5: Values (3\ and (5 2 - 

Now we try approach based on similarity of nucleation conditions [Sj. We 
start with the infinite upper boundary in formulas for f3\ and /?2- 

At first we shall show the result for renormalization of dispersion when 
for a we take the result from the two cycle model. The result is shown in 
figure 9. 

At p = 0.64 the value of ip will be 0.46, at p = 0.68 the value of ip will be 
0.43. 

When we put a = 0.92 which is the result in the precise solutuion we 
came to the result shown in figure 10. Practically nothing has been changed. 

Now we shall restrict the upper limit of integration by z = 1 (this is the 
current moment of observation). At first we draw the number of droplets as 
function of p in figure 11. 

We see that the number of droplets has minimum at p = 0.43. This value 
has to be chosen as p in calculation of dispersion. 

Then the value of relative square of dispersion can be calculated and it 
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Figure 6: Ratio fii/fo 



is drawn in figure 12. 

We see that at p — 0.43 the value of dispersion is near 0.64 and it is 
far from the expected value. It resembles the calculation without similarity 
taken into account. 

It is interesting to note that in diffusion regime the values of relative 
dispersion ip calculated in the model with the upper boundary 00 and with 
the upper boundary 1 are practically the same. This is true for moderate p 
(more concrete p < 0.7). 

An explanation of this property is the following: under the free molecular 
regime of droplets growth the droplets appeared in the first moments of time 
grow so rapidly that it was very important to know until what sizes (or until 
what moment) they will grow. So, the duration of the nucleation period was 
very important. Under the diffusion regime of growth this effect is much 
more weak. So, the position of the boundary isn't so essential. 

For p < 0.5 the value 1 — (3 /a is also close to these values. 
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Figure 7: The relative square of dispersion in the two cycle model. 

The explanation lies also in collective character of stochastic effects ini- 
tiated under the diffusion regime of droplets growth. 

The question to decide is whether we have to take p from extremum for 
a calculated with the upper boundary oo or with the upper boundary 1. 
From the first view it seems that the boundary 1 is preferable because f3i 
and p 2 are calculated with the upper boundary 1. But the extremum for a 
corresponds to the value a = 0.815. This value is smaller than a calculated 
in the first iteration, i.e. a\ = 0.8773 and smaller than this value calculated 
in the precise solution a pr = 0.9292. So, at extremum (it is a minimum) of 
the model with the upper boundary 1 the deviation from the precise value 
attains maximum. At the upper boundary oo the value of a is always bigger 
than 1 which is greater than ol\ and a pr . So, here the extremum is really the 
closest value to precise solution. Then it has a real physical meaning. So, it 
is reasonable to take p corresponding to the extremum of a with the upper 
boundary oo. 
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Figure 8: The relative square of dispersion in the two cycle model. Account 
of the change of the droplets number. 

Calculations give p = 0.65 and ip = 0.49. This value is calculated on the 
base of f3\ and $2 with the upper boundary 1. The value ip = 0.49 is rather 
close to result of numerical simulation ip = 0.44. Certainly, the relative error 
is one tenth and the methods based on the explicit two cycle models cannot 
give more precise results. 

One can note the following curious result. Evidently when we calculate 
ip for p = 0.65 with the upper boundary 1 we come to ip = 0.46 which is 
practically the necessary result. But this coincidence is no more than an 
occasional one. When we try the same procedure for the free molecular 
regime we come to ip = 0.55 (here ip = 0.55) which is far from numerical 
result ip = 0.64. 

Having estimated the value of deviation of ip = 0.49 from ip = 0.44 one 
has to note that it is much smaller than the deviation of ip — 0.59 of the 
value 1 — (3 /a in its minimum from ip = 0.44. So, the modification made 
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Figure 9: The relative square of dispersion in the similarity model. Account 



of the change of the droplets number in two period model, 
here really have a real sense. 

3 The effect of growing volumes of interac- 
tion 

What shall we do to take into account the growth of the volume interaction 
in diffusion regime of growth? It is necessary to rewrite formula (j3J) • We can 
substitute in exp(— ^"^^ ) the denominator 2iV\ by another denominator 
which corresponds to the growing volume of interaction. Then Ni has to be 
substituted by NiConst(j — i) 5 ^ 2 with some known const. It means that the 

characteristic half-width \[W\ has to be substituted by <J 1 Niconst(j — i) 5 / 2 . 
One can not apply directly this substitution in exp(— ) but only ap- 
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Figure 10: The relative square of dispersion in the similarity model. The 
droplets number is taken from precise solution. 



proximately in expression for a\ which now will be 

with some known const and an appropriate renormalization. 
In diffusion regime 



and there is absolutely no effect. It means that kinetics of nucleation will be 
absolutely different. It is described in jlj. The most strong effect will take 
place in the free molecular regime where j^- + 1 — 5/2 = 3/2 

Now we shall estimate the effect of the dispersion variation. We see 

D (oo) _ ^(oo)^ _ _( 2 (^_ + 1) / dx dz{z - x)^ +1 - 5/2 exp(-^ 5/2 ) 
a 1 — s Jo Jp 
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Figure 11: The number of droplets until the current moment. 




( 1 +1) 2 I" 'dxiT 'dz(z x)T^ +1 - 5 / 2 exp( z 5 ' 2 )) 2 )) 

1 — S JO Jp 



Note that the last relation has to applied to the subsystem with the size 
\fADtn where t n is the duration of the nucleation period. Then we have to 
take the superposition of these independent subsystems. 

We see that in the majority situations we can approximately neglect (j — 
i)T^ +1 ~ 5/2 w (j - i)° = 1 and see that then 

£,(oo) ^ jy-(oo) 

and there is no effect. 

Nevertheless the deviation in the number of droplets from the value cal- 
culated on a base of averaged characteristics exists. We have to divide the 
whole system into subsystems of the size y/ADt n , then determine the number 
of droplets in this subsystem N^ s obtained on the averaged characteristics, 
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Figure 12: Relative square of dispersion up to the current moment. Account 
of the change of the droplets number. 

then determine the deviation 5N^ t bs from the values calculated in TAC and 
then the whole deviation in the whole number of droplets N to t from the value 
based on the averaged value N tat will be 



These procedures completely solve the problem. 

One has also to note that the use of monodisperce approximation ,2] leads 
to the absence of the problem of growing volume because the characteristics 
governing the nucleation period have been determined only in one moment 
of time. The problem here is a justification of monodisperse approximation 
to the calculation of fluctuation effects. 
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4 Appendix 



Here we shall present the two cycle model with a fixed boundary in kinetics 
of the metastable phase decay under the free molecular regime of the droplets 
growth. This question was completely investigated in jS] but during the last 
two years our view on this problem was slightly modified and the modern 
state is described below. 

Systematic investigations of the first order phase transitions are per- 
formed since Wilson The classical theory of nucleation |8 a gave expres- 
sions for all main characteristics of stationary process of nucleation. This 
allowed to investigate a global picture of the phase transition. A set of pa- 
pers was devoted to model pictures of the global kinetics of nucleation. 
Here we shall also consider the global picture of the phase transition. One 
can note that all cited publications were based on averaged nucleation 
rates. Here we shall consider stochastic appearance of embryos and outline 
the stochastic manner of appearance. 

Recall briefly the main features of phase transition. Suppose that in 
initial moment of time there exists a metastable state. Then the embryos of 
a new phase begin to appear in the metastable system. The average rate of 
appearance is given by e 8j. Then the embryos begin to grow and to consume 
the vapor, metastability falls and the rate of nucleation, i.e. the rate of 
appearance of new embryos falls also. The vapor consumption occurs in a 
time scale in a very sharp avalanche manner. 

It is clear that during the nucleation period the new supercritical forma- 
tions of the new phase appear with some fixed probability, but they appear 
in stochastic manner. So, the stochastic appearance of relatively big number 
of droplets leads to very rapid stochastic consumption of vapor. Stochastic 
appearance of relatively small number of droplets leads to delay of the vapor 
exhaustion and to excess of droplets appearance in next moments of time. 
But it seems that this excess can not compensate the opposite effect of the 
absence of embryos in first moments of time. It seems that the total number 
of droplets will differ from the average value. This is the naive reason why 
stochastic effects of nucleation have to be taken into account. 

In a system with macroscopic sizes due to a giant value of the Avogadro 
number there appears some rather big number of droplets. It allows to use 
the averaged characteristics to construct kinetics of a nucleation process. 
Precise kinetic approach based on averaged characteristics is described in 
In |3| the time evolution is completely constructed. 
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After the formulation of integral equations (see |3J) one can introduce 
"elementary intervals of nucleation" - the intervals where the state of the 
system changes negligibly small. In macroscopic systems the total number 
of droplets is so big that at every elementary interval there appears a great 
number of droplets AN. On the base of traditional thermodynamics one can 
state that the relative fluctuation 5 AN /AN of droplets formed at elementary 
interval is small and has an order of (AiV) -1 / 2 . This remark completely 
solves a problem of justification of nucleation description based on averaged 
characteristics. 

Since it is possible to extract elementary intervals where thermodynamic 
parameters and the nucleation rate have small variations there is no need to 
take care about the stochastic corrections. 

In experimental investigations one can not study quite macroscopic sys- 
tems because the most popular data is the number of droplets and to fulfill 
the calculation of this number one can not have too many droplets 3 . After 
the theory based on the averaged characteristics has been presented it be- 
came possible to investigate the stochastic effects in kinetics of nucleation, 
i.e. the effects of stochastic appearance of droplets. Recently there appears 
a set of papers jTU], PQ, [H] where a stochastic effects (the effects of fluctua- 
tions of droplets formation) were described and investigated. One can extract 
two approaches which were formulated in and in pQ. Although pQ, [TT| 
were written by the same authors of [TO] these authors didn't hesitate that 
the theories formulated in [JJ , [TU] gave different results. So, at least one 
has to analyze approaches [TU], [JJ and decide whether there is a true result 
among these approaches and how it can be used to construct the adequate 
description of the nucleation process. 

One has to specify a formulation of the problem of our investigation. 

In P], ^UJ it was proposed to establish corrections to the total number of 
droplets iV appeared in the system. It was supposed that these corrections 
are functions of N. To demonstrate the error of this approach it is sufficient 
on one hand to take two identical systems then to calculate them separately 
and to add results or on the other hand to calculate correction directly for 
the total system. The results will be different. 

One has to determine a real volume to which one has to refer the number 
of droplets. It is simple to do with the help of results from [1J. In that 

3 One can not simply calculate the infinite number of droplets. The upper limit of the 
number of droplets which can be noticed on a photo image is about several dozens. 
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paper the kinetics of nucleation for diffusion regime of droplets growth was 
constructed. It was shown that a solitary droplet perturbs vapor up to dis- 
tances of an order y/4Dt, where D is a diffusion coefficient, t is a time from 
a moment of droplet formation up to a current moment. As an estimate one 
can take as t a time t 1 of the whole nucleation period duration 4 . 

The last remark allows to give a new definition of the volume V e i where 
the number N e i of interacting (mainly through the vapor exhaustion) droplets 
is formed. Namely this value has to be regarded as a volume of a system in 
approach [TU], [Tj. This volume is 

V el = 4vr(4 J Dt 1 ) 3 / 2 /3 

If the sizes of the system are smaller than this value one has to take the 
volume of the system as this value. But such a small system can be hardly 
regarded as a macroscopic one. At least one has to analyze carefully the 
boundary conditions for the system. 

Naturally, the droplets appeared at different moments of time perturb 
initial phase up to different distances. So, one can regard above formulas 
only as estimates. Some more rigorous equations can be found in [13j. 

The number of droplets N e i isn't too big as N is. So, an analysis of 
stochastic effects has a real sense. It is interesting now to get all correction 
terms which are ascending with the number of droplets (but not only a leading 
term). To solve this task one has to modify approaches from pQ, |10j . 

Complexity of this problem appears here because one can not directly 
use equations based on the theory with averaged characteristics. In pQ, JUJ 
some properties of solution of the theory on the base of averaged characteris- 
tics (TAC) were the starting points for constructions. This supposition was 
adopted without any justification. So, at first we have to decide whether it 
is possible to start with TAC. 

We shall consider the situation of decay of metastable phase. The new 
dimensionless parameter - the number of droplets destroys the universality 
observed in [3] for the theory based on averaged characteristics. Moreover, 
it is difficult even to formulate the system of equations. This radically com- 
plicates the problem. 

The possibility to use the effective monodisperse approximation formu- 
lated in [2] was used in [TU] without any justifications. Generally speaking 
this property can not be directly used to calculate stochastic corrections. 

4 The nucleation period is a period of intensive formation of droplets. 
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One has to analyze whether this conclusion leads to an error made in [TU] . 
Here we have to formulate more correct constructions. 

Both approaches from pQ and from ^U] declared that they used the fol- 
lowing property observed in [3] in frames of TAC: 

"The droplets formed at the beginning (i.e. at the first half) of the nucle- 
ation period are the main consumers of vapor" . 

This property is valid but it is substituted in reality in [TU], PP by the 
following statement: 

The main source of stochastic effects are the free fluctuations of droplets 
formed at the beginning (i.e. at the first half) of the nucleation period. They 
govern the fluctuations of all other droplets. 

The last statement seriously differs from the first one. To get credible 
results it was necessary to balance the fluctuation effects from the first half 
of spectrum by the corresponding reaction of remaining part of spectrum. 
Then, at least this approach needs some special justifications. So, one has 
to use some new constructions which are presented below. 

The application of the model approximation which was in reality done 
in PP, ^JJ will lead to some errors. But due to universality of solution j3] 
the errors can not be be too big. Qualitatively everything is suitable, but 
precision will be not so high. 

The same conclusion will be valid for any approach based on some model 
behavior of supersaturation (justification is valid for a vapor consumption 
in TAC, but not for stochastic effects). Namely, in pQ was used an artificial 
approximation where at the first half of the nucleation period the conditions 
of nucleation are the ideal ones and at the second half the conditions of 
nucleation are governed by the droplets formed at the first half. Here in 
current paper the final result will be more precise and it will be not based on 
rather spontaneous artificial choice of some parameter as it was done in pQ 
where this parameter was put to to 1/2. In [1] it is supposed that until some 
moment (it is chosen in pP as a half of nucleation period) the droplets are 
formed under ideal conditions and namely these droplets determine a vapor 
consumption. In reality this approach taken from [T2] (page 310) was used 
in pP in slightly another sense. It is supposed that droplets formed during a 
first half of nucleation period are the main source of stochastic effects. The 
last statement was not justified in pP and it is rather approximate. The 
relative correctness of a result was attained due to specific compensation of 
different errors of approximations used in pp. It is necessary to stress that 
the the mentioned model was used in ^2] to justify a strong inequality and 
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the high precision of constructions was not essential. But in PQ this model 
was the source of numerical result. 

All arguments listed above lead to necessity of reconsideration which will 
be made in this publication. A plan of narration will be the following one 

• Having considered the interaction of stochastic deviations of the num- 
ber of droplets appeared during the elementary intervals of nucleation 
we shall see that stochastic effects are at least moderate. 

• A moderate scale of stochastic effects allows to seek the solution on 
the base of the theory with averaged characteristics. But we have to 
take stochastic effects from all droplets formed during the nucleation 
period. 

• The possibility to take into account the influence of stochastic devia- 
tions of all droplets can be provided by the property of the self similarity 
of nucleation conditions during the nucleation period. This property 
can be considered in two senses - 1) in the local differential sense and 
2) in the integral sense in frames of the first iteration in the iteration 
procedure in TAC 15] . The local property will be used in justifica- 
tion of the smallness of stochastic effects and the integral property will 
be used to calculate the stochastic corrections of the whole nucleation 
period. 

All analytical results will be checked by computer simulation and a coin- 
cidence will be shown. 

All mentioned constructions will be valid for an arbitrary first order phase 
transition. The law of droplets growth here will be a free molecular one, then 
the linear size of droplet grows with velocity independent from its value. 
Consideration of other regimes can be formally attained in frames of the 
current approach by some trivial substitutions, but one has to take into 
account that the new regime requires new approaches to construct nucleation 
kinetics as it is shown in 4J. So, one can not agree with JU] where it is stated 
that one of results is an account of stochastic effects in a diffusion regime 
of droplets growth. This effect has to be taken into account principally in 
another manner by application of methods presented in [Tfj . 
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4.1 Some characteristic features of decay kinetics 



We begin with the theory based on averaged characteristics. It is supposed 
to be known [3] , that the supersaturation £ behavior can be determined after 
certain renormalizations by the following equation 

ip(z) = / dx(z — x) 3 exp(— ip) 
Jo 

on function ip which is the relative renormalized deviation of supersaturation 
from the initial value. Variables x and z can be considered as equivalent 
ones. A good approximation for solution and for a distribution of the droplets 
number over time (or over liner sizes of droplets) which is proportional to 

exp(-0) is 

/i = exp(-**/4). (4) 
The form of f\ is given by fig. 13. It is seen that approximately at zq = 1.25 
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Figure 13: A form of approximation for size spectrum 
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the nucleation period stops. 

Approximation (J3J) has rather high precision It is based on the fol- 
lowing law of substance accumulation 



pz pz 

G = z 4 /A = / (z - xfdx = / g(x)dx, g = (z - xf = p 3 
Jo Jo 

for the renormalized number G of molecules in a new phase. For any moment 
t or z a function g has one and the same form. We shall call this property as 
a "similarity of nucleation conditions". We see that every time the droplets 
formed at the last third of a period from beginning of nucleation until a 
current moment will accumulate a negligible quantity of substance. The 
relative quantity of the substance there has an order of (~ 1/27) and is so 
small that even if there will be fluctuations the quantity will remain small. 

From the form of f\ it is seen that until z= = 0.7 zq all droplets deplete 
vapor rather weak. It will be important for future analysis. 

Another important property is the possibility to describe kinetics in frames 
of TAC with the help of monodisperce approximation (see [2 ). The men- 
tioned property of g allows to use a monodisperse approximation j2] not only 
at the end of nucleation but in every moment of the nucleation period j2] • Let 
t(G) be the moment when there are G molecules in droplets (in appropriate 
units). An application of the monodisperse approximation [2] leads to 

G « N m (z)z 3 

where N m (z) is the number of droplets born until z/A 

N m (z) = N(z/A) « z/A 

and 




for z which are essentially less than those corresponding the end of the nu- 
cleation period (in reality z < z=. 



4.2 Interaction of arbitrary stochastic fluctuations 

The account of fluctuation interactions at every moment of time is important 
in justification of the smallness of stochastic corrections. Generally speaking, 
one can not take consideration only at the end of nucleation period without 
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justifications that the characteristic features of nucleation kinetics remains 
the same at every moment of time. 

The arbitrary value of z in the period of nucleation corresponds to the 
arbitrary value of the number of droplets in a liquid phase G (in renormal- 
ized units which will be used (see the value G will belong to interval 
(0, 1)). Now it will be possible to repeat in some features the approach from 
[TUj with arbitrary parameter G instead of 1 (in renormalized units, before 
renormalization it would be 1/T (see jS]))- The sense of a difference from ^U] 
is consideration of an arbitrary moment z instead of the end of nucleation. It 
is very important because allows to take into account all fluctuation effects 
during the nucleation period. 

Let t(z) be the current moment of time (z is the coordinate of the spec- 
trum front, actually t is proportional z). We suppose that before az (a is 
some parameter) droplets are formed without mutual influence and one can 
write Poisson's distribution. This is the first group of droplets. The sec- 
ond group of droplets are all other droplets formed until the time moment 
z. A natural restriction on a appeared, namely a < 0.7. We shall take 
a > 1 — 0.7 = 0.3 also for the purposes explained below. We suppose that 
the influence of other droplets on its own formation is negligible (this fol- 
lows from 2 * 0.7 = 1.4 > 1.25 and from notation made above about the last 
third of nucleation period). Then it is possible to write Poisson's distribution 
for the second group of droplets, but with parameters depended on stochas- 
tic values - characteristics of the droplets distribution from the first group. 
Rigorously speaking one has to use the first four moments of the droplets 
distribution in accordance with [3 , but for simplicity we shall use here only 
the zero momentum. As a compensation for this simplicity we has to use 
here only a = 1/4 which corresponds to the applicability of monodisperse 
approximation. But due to the arbitrary value of G this is quite sufficient 
for our goal. 

Certainly, one can not state that precisely the first fourth of the spectrum 
plays the main role in vapor consumption. So, we have to consider interac- 
tions of stochastic fluctuations for all a which aren't too small in comparison 
with 1/4 as well as 1 — a isn't too small in comparison with the same 1/4. 
This will be done below. 

At the next step of consideration one has to come from Poisson distribu- 
tions to Gauss distributions and integrate them with account of connection 
between stochastic parameters of embryos formation from the first group and 
parameters of distribution from the second group. The same was done in ^U] 
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but only for a leading term. 

Why it is necessary to get all ascending terms in corrections? The ex- 
istence of at least one coefficient with a big absolute value means (see the 
theory of Chebyshev polynomials) that there exist a size of a group where the 
interaction will lead to the big effect. Then it would lead that the interaction 
in these groups will be the real driving force of the process and these groups 
can be regarded as quasiparticles. 

Contrary to ^U] we shall take into account all correction terms which 
come from transition from Poisson's distributions to Gauss distributions and 
corrections for nonlinear connection between the group distributions. We 
shall take all terms which are growing when the total number of droplets 
grows. 

We get the following result for droplets distribution 

P = P 0O (l + y) 

where 

°°" l 27riV(3a + l) j Pl 3a + 1 2 j 
N-N 



N 



N - some stochastic value of the total number of droplets, iV - the mean 
value of droplets and y is the correction for spectrum. 
At =1/4 we get 

^^ D ^ 087D -- W2e9 > + i -l + m- 6 D ' 

85903 ^ 4 65399569 ^ 6 2 
~~ 12446784 s + 10978063488 s ^ 

where 

s = i/Vn 

is a small parameter of decomposition. To get all ascending (with N) cor- 
rections we must fulfill decomposition until s 2 . 

Why it was necessary to get all ascending terms in decomposition? The 
answer lies in specific sequential influence which can be observed in nucleation 
period. The droplet appeared in the first moment of time forms condition for 
the embryos appearance in the second moment, then the embryos appeared 
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in the first and in the second moments form conditions for further appear- 
ance, etc. So, if there would be some N at which corrections will be big 
then immediately namely this value will be the crucial value for all kinetics. 
Fortunately here there is no such effect and, thus, we can take the theory 
with the averaged characteristics as the base for further constructions. 

At arbitrary a we get for P/Pqo (here is the limit at N = oo) the 
following expression 

P/Poo = l + W 1 S + W 2 S 2 



Here 



Wl = _I((486/ 12 + 486P 2 / 12 - 972/ n P 2 + 324/ 10 - 648P 2 / 10 + 756Z 9 



+810/ 9 P 2 - 459/ 8 + 297P 2 / 8 - 135/ 7 P S 2 - 
387/ 7 +90/ 6 P 2 -153/ 6 +27/ 5 P 2 -213/ 5 -66/ 4 P 2 -3/ 4 +27/ 3 -3/ 3 P 2 +16/ 2 P 2 -3/ 2 + 
9/ + ID 2 + D 2 )D S )/((1 + 3Z 2 ) 3 (Z + 1)(-1 + I 2 )) 



where 
and 



1 = \[a 

w 2 = w 02 /((l + 1) 2 (-1 + Z 2 )(l + 3Z 2 ) 6 / 2 ) 

3 

9o = -i - i( + mr _ ™e _ m * _ «Lp + + *v 

yo 12 6 4 2 2 2 2 

153 z io _ 39 /5 _ 67 /4 _ 225 /9 _ 3 2 _ 81 u 243 13 _ 243 /15 _ 17 
2 26222 2 2 6 

q 2 = + ^Z 6 + 117/ 10 + S 11 + 

+ ^ 5 + ^ 14 _ 15309 zl9 + 2079 p 
2 2 2 

+31Z 8 - 8748/ 18 + 42Z 9 + ^> + — I 7 + l^Zp + 6561 p 

2 2 8 

+ !^IZf _ 3159/ i 6 _ 8262/ ir 
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„ = - ™(» + 7047/- - IHp 

+ 3M i8+ 105 ;11 _104|47 i22 _3 i4 + 170Hi20 

_ 5103 p + _ + 7%5i " -m»- ™l" + SW» - + ^« 

6 4 4 3 6 4 

1 ;3 44469, 24 1 , 2 6903 , 16 74 ;9 

H Z 3 H Z 24 Z 2 Z 16 H Z 9 

12 4 12 4 3 

40419,0 17253 l9n 24057, 99 ~, 6561 l9fi 

g 6 = -13122/ 25 Z 18 + Z 20 I 22 + 4374Z 24 + 1 26 - 

8 2 2 2 

48843 21 _ 623 /8 59 /6 _ 3609 14 

2 24 36 4 

_57 411 11 ± _ 29 _ 345 1 

8 4 24 72 12 4 36 

13 )5 65,„ 2421 , 13 7857, 15 

— Z 5 H Z 9 H Z 13 Z 15 + 

36 4 4 4 

— Z 17 + 28431Z 23 + 7047Z 19 + — Z 16 
4 8 

Having integrated this expression we get corrections to droplets number. The 

term at s gives zero after integration and the first non zero correction has 

an order of ~ s 2 and doesn't depend on the total number of droplets. A 

coefficient at s 2 has at a — 1/4 a value 

d = 311/3024 < 1 

At arbitrary a a coefficient at s 2 in correction for the total number of droplets 
will be 

d a = —[108a 6 + 540a 11/2 - 72a 5 - 930a 9/2 - 336a 4 + 

713a 7/2 + 158a 3 + 4a 2 - 203a 5/2 - 6a + 39a 3/2 - 3a 1/2 ]/[a 3/2 (l - a) 2 (l + 3a)] 

It will be interesting to compare results with and without corrections from 
transition from Poisson's distribution to Gauss distribution. So, we consider 
now this case. At the leading term there will be no change. At correction 
terms we have 

y = — DJ-2331 + 289D 2 )s + (—D 2 + 

y 74088 V sJ V 196 s 
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732037 dA 24137569 ^ 6 13 2 
12446784 s + 10978063488 s ~ 3fT 
d = -37/126 

d a = -l(648a 11/2 - 216a 5 - 1062a 9/2 + 108a 4 + 753a 7/2 - 

30a 3 - 195a 5/2 - 12a 2 + 19a 3/2 + 6a - 7a 1/2 )/(a 3/2 (l + 3a) (1 - a) 2 ) 

It is seen that these corrections are small. At arbitrary a except too small 
ones and those close to 1 (these values are unreal) we get values shown at 
fig. 14 
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Figure 14: Dependence of corrections, i.e. of d a on a 



A point curve shows corrections with transition from Poisson's distri- 
bution to Gauss distribution taken into account, a line shows corrections 
without transition from Poisson's distribution to Gauss distribution. Both 
corrections have one and the same order and they are small. 
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We see the plateau for all a except a C 1/4 and I — a <C 1/4. So, here the 
smallness of corrections due to the interactions of stochastic peaks is quite 
clear. 

But still the small value of corrections can lead through sequential in- 
fluence to essential change of the total droplets number. Conclusion about 
the smallness of correction in the total number of droplets can be made on 
the base of specific kinetic features of the process (see "the effect of the first 
droplet") which will be done later. 

One has to note that we have not taken into account corrections from 
transition from summation to integration. This is definitely required by 
discrete character of droplets number. It is made due to reasons formulated 
below. Really, we have at transition from summation in formula 

P(N) = Y,P 1 (N 1 ,N 1 )P 2 (N - N, - N^N^NJ) 

Ni 

to integration 5 

P(N) = J dN^A, N 1 )P 2 (N — N 1 — N 2 , A^i)) 

to use the Euler-McLorrain's decomposition. It brings to asymptotic serial, 
which can be included into a final answer. 

This is the formal solution of the problem. But discrete character in 
nucleation isn't so trivial. The process of vapor consumption can not begin 
without the first droplet. The system will wait for droplet as long as it 
will be necessary. It shows that discrete effects are complicate and require 
a separate publication. At least one has to put the initial moment at the 
moment of the first droplet appearance and then to consider condensation 
with the substance in the first droplet calculated explicitly as ~ z 3 

This property will be called as "the effect of the first droplet". Here it 
leads only to the small effects, but in the situation of smooth variation of 
external conditions the effect can be greater. 

5 Here N is the total number of droplets, Ni is the stochastic number of droplets in the 
first group, N2 is the stochastic number of droplets in the second group, N\ is the mean 
number of droplets in the first group, N2 is the mean number of droplets in the first group, 
which is a function of Ni. The value N remains stochastic number. 
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4.3 Self similarity of Gaussians 

To use Poisson's distribution for the first group of droplets one has to make 
the following notation. Really nucleation conditions for the first group don't 
differ from the whole group. So, for distribution for the first group one has 
to take distribution Pi with reduced half-width. But one can not attribute a 
half-width to Poisson's distribution. That's why we considered effects with 
and without corrections from transition from Poisson's to Gauss distribution. 
So, we can use Gauss distributions as initial ones. For Gauss distribution 
one can easily reconsider the half-width. Then for Pi one can take 

Pi ~ exp(-(X\ - N 1 ) 2 /(2bN 1 )) 

where b is a renormalization coefficient. Distribution P2 remains previous 

P 2 ~ exp(-(A> 2 - N 2 ) 2 /{2N 2 )) 

where N 2 is given by 

„ 1 2 2 14 3 35 4 91 5 

N 2 = (1 - -s + -s 2 s 3 + s 4 s 5 - a)N 

y 3 9 81 243 729 ; 

where N is a mean total number of droplets, 

Ni-aN 
S ~ aN 

is a small parameter of an order N~ 1 / 2 
After integration one comes to 

GXp(_ 2(9a + b - 6ba + 96a 2 - 9a 2 ) } 

where 

d= (N-N)fVN 

The half-width of the distribution Pi must be equal to the half-width of 
P, which leads to 

b = 9 f~ a) (5) 
-9a 2 + 15a - 1 K ' 

This value is drawn in fig. 15. 
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Figure 15: Function b(a). 

Now we shall consider effects of renormalization. The ratio of corrections 
with renormalization and without renormalizations is given by 

_ 1 V-9a 2 + 15a - ly/3a + 1 
9 ay/1 — a 

and it is shown at fig. 16 

For all reasonable values of a the last ratio is approximately 1. At a = 1/4 
we get 1.0041 Here the effect of similarity of nucleation conditions doesn't 
lead to remarkable effects in the change of half- width. But here the smallness 
of change takes place only because of application of monodisperse approxi- 
mation and later this change will be essential. 

Instead of taking into account all moments of distribution we can directly 
calculate the effects on the base of explicit form of spectrum in frames of 
iteration procedure. 

The result of the previous consideration is the following: we have proved 



41 



1.6- 
1.4 : 
1.2 : 



0.8- 
0.6 : 
0.4 : 
0.2 : 



0.2 0.4 0.6 0.8 

a 

Figure 16: Ratio of halfwidths 

that stochastic effects are small. Beside this we have demonstrated how to use 
the similarity of nucleation conditions. Here it doesn't lead to any remarkable 
effect, but later under the smooth external conditions this method will lead 
to some essential numerical corrections. 

One has to stress that in pQ a linearization over deviation of the droplets 
number from the mean value was made. It was necessary to perform calcu- 
lations. But in the linear theory one can not get the deviation of the mean 
value of result from the value calculated on the base of averaged character- 
istics. So, the attempt to get deviations in the mean value of droplets in £Q 
is senseless. 

4.4 Calculation of dispersion 

Since the mentioned combination of gaussians is also a gaussian characterized 
by a mean value and by a dispersion one has to determine these values for 
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the distribution of the total number of droplets. As it was shown above the 
mean value of the droplets number is practically the same as it is prescribed 
by TAC and we are interested in dispersion. Now we shall calculate this 
value properly. 

The most advanced approach to solve this problem 6 was suggested in [I] . 
But even this approach has many disadvantages and we need to reconsider 
it. 

We shall characterize a droplet by a linear size p which is the cubic root 
of its molecules number. Its velocity of growth at fixed supersaturation does 
not depend on p. 

Decomposition of a whole interval of nucleation into elementary intervals 
is connected with some difficulties. An elementary length A according to pQ 
must satisfy two requirements: 

1. A number of droplets formed during elementary length must be great. 

2. An amplitude of a spectrum has to be approximately constant during 
an elementary interval. 

It is clear that the second requirement can not be satisfied. Stochastic 
deviations of an amplitude leads to the violation of the second requirement. 

We shall apply the second requirement not to the stochastic amplitude 
as it was stated in pQ, but to the averaged amplitude. Then the second 
requirement is: 

An averaged amplitude of a spectrum has to be approximately constant 
during an elementary interval. 

Now we shall see the evident illegal consequences of the approach from 
PQ. "Stochastic" amplitudes fi are introduced in pQ as 

f. = ^ 

A 

where iVj is the number of droplets formed during A. It isn't the height of 
spectrum but simply the renormalized value of the droplets number appeared 
during this interval. An expression for the number of molecules in droplets 
formed during interval number i at a moment tk (or Zk) (it means that now 
we are at interval number k) with approximately constant rate of nucleation 

6 Certainly, already expression © gives the decrease of the Gaussian halfwidth and can 
lead directly to dispersion where b is smaller than a standard one. But b essentially depends 
on a. For 0.2 < a < 0.9 one can see the approximate formula b = 0.6 — (a — 0.5) ^gE§l 
and the absence of concrete value for a doesn't allow to get a concrete value for b. 
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fi would be the following one 

/ Sip dp = -fi(x k _ i+1 - x k _i) 

Jx k _i 4 

Namely, this equation was derived in £Q and forms the base of further con- 
sideration. Here Xk-i is the coordinate p of the droplet which was born at Z{ 
at a moment z k (it corresponds to the definition x = z — p) The difference 
between forth powers corresponds to a constant amplitude of spectrum. It 
is wrong and then eq. (12) in pQ and all further equations are not correct. 

But the is no necessity to use such a way to make an account of the 
number of molecules in a new phase. It is absolutely sufficient to take the 
following expression 

/ fip'dp = Nixt^ « N iX s k _ i+1 

which is valid at k — % 3> 1. In a whole quantity of substance it is sufficient 
to take into account only droplets with k — i ^> 1. The relative weight of 
dismissed terms will be small. 

Then for the total number of molecules in droplets at interval number k 
we have the following expression 

Q k = ]T 

i=i 

where corresponding coordinate. This expression can be rewritten 

as 

k 

Qk-Y^N^ik-if 
i=i 

This representation is important because now the note in [I] after eq. 
(15) isn't necessary. That note stated that the probability for iVj to deviate 
from the number iVj of droplets calculated in TAC under the supersatura- 
tion formed by stochastically appeared droplets in previous intervals is very 
low. That note is doubtful because namely these deviations are the base for 
stochastic effects. Now there is no need in this note. 

The next step is to build [lj a two cycle construction for nucleation period. 
During the first cycle the main consumers of vapor appeared in a system 
and during the second cycle they govern a process of formation of all other 
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N tot = - — ,\ = 1.2818 



droplets. In [T] it is supposed that during the first cycle a vapor depletion is 
negligible and during the second cycle new droplets are absolutely governed 
by droplets from the first cycle. Now we shall analyze an effectiveness of such 
procedure. 

In TAC the corresponding evolution equation will be 

ip(z) = / (z — x) 3 exp(— ijj(x))dx 
Jo 

The first iteration [3] is practically a precise solution and it gives the number 
of droplets 

14^W2 
4 T(3/4) 

A model solution requires that until z = p there will be no depletion of 
vapor and then only the droplets formed before z = p will consume vapor. 
Then for a total number of droplets we have an expression 

/oo 11 
exp(--x 4 + -(a; - pf)dx 

A ratio q = N tot apw /N tot is given in fig. 17 

Always N tot appr is greater than N tot . The value of minimum corresponds 
to p ~ 0.78 which is 55 percent of the total length of spectrum. We can 
stress the smooth dependence N tot appr on p. 

It is clear that in pQ the value of parameter of separation into two cycles 
was not chosen in a good style (at least from the point of view of TAC). It 
corresponds to p — 0.64. 

Now we shall study the probability Pk of formation of stochastic number 
iVfc of droplets at the first k elementary intervals. Our constructions now 
resemble |IJ but there is one essential difference. We have no necessity to 
linearize expression with respect to (Ni — N^/Ni, where iVj is a stochastic 
number of droplets formed at interval i, Ni is a mean number of droplets 
formed at interval % (it is a function of stochastic numbers of droplets at 
preceding intervals). This linearization can not take place because a ratio 
(Ni — Ni)/Ni can be zero or can attain huge value (with a low probability). 
It is more simple and more justified to linearize expression on J2iPi(N — 
N)/Ni where pi is a linear size of droplets formed at interval i (all of them 
have approximately the same size). Really, due to summation the relative 
variations of J2iPi(Ni — Ni)/Ni are much smaller than variations of (A^ — 
Ni)/Ni. 
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Figure 17: The ratio of mean numbers of droplets 

Variations of (iV, — Ni)/Ni would be small only at very big numbers of 
droplets N tot . One can get 

Ntot ~ MNi 

M is a number of elementary intervals. So, the theory with linearization 
proposed in pQ would be well justified only in a region where the result can 
be obtained on the base of averaged characteristics. The internal contra- 
diction between the big number M and the smallness of fluctuations in the 
elementary interval in pQ is obvious. 

The linearization proposed here is much more weak than in p. But it 
leads to the analogous numerical expressions as in pQ. So, restrictions from 
PQ are not necessary. 

For dispersion of the total distribution the result proposed in pQ was the 
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following 

D°° = N°°(l - -) 
a 

where 

P = Pl-02 

poo 

A = 8/ dU rfr(e-r) 3 exp(-e 4 ) 

Jl/2 JO 
poo poo /"1/2 

/3 2 = 16 / d£ / d77 / dr(r - 3 (r - V ) 3 exp(-^) exp( V) 

Jl/2 Jl/2 JO 

a = dx exp(—x 4 ) 
Jo 

In the two-cycles construction the value of a, which is proportional to the 
total number of droplets has to be reconsidered and recalculated on the base 
of two cycles. Then we have to use instead of previous a = «o a new value 

poo 

a = oil = 1/2 + / exp(-x 4 + (x - l/2) 4 )cfa 

A/2 

In our approach we shall use parameter k of separation 7 of two cycles and 
we shall calculate a± as 

POO 

«i = k + / exp(— x 4 + (x — k) 4 )dx 

Jk 

Then according to fig. 17 we see that the ratio ao/ ot\ is greater than 1 and a,\ 
is greater than a - Here we see that two-cycles construction is approximate 
one. Then result for 73°° will differ from the number published in [I] and will 
be (here one has to put k — 1/2) 

7300 = jyooQ.69 

instead of 

£,°o = ^ooQ.67 

as it is written in pQ. 

Numerical simulations show that the value is one tenth more than a 
real result. So, the new theory is necessary. 

7 Here for simplicity we use k instead of p/4 1 / 4 . 



47 



Now it is necessary to decide what k shall we choose. At arbitrary k the 
expression for (5 will be the same but for f3\ f3 2 we have 

poo rk 

A = 8/ dd rfr(e-r) 3 exp(-e 4 ) 

Jk JO 
poo poo pk 

(3 2 = 16 / dd dr, dr(r - £) 3 (r - r,f exp(-£ 4 ) exp(-r ? 4 ) 

Jfc JA; JO 

We have also to reconsider expression for a. 

After calculations we have for dispersion as function of k the following 
fig. 18 
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Figure 18: Relative deviation in dispersion as function of k 

A minimal dispersion will be at k = 0.6. This value approximately equals 
to k — 0.55 which ensures minimum of ol\ where the result is the most close 
to the real value in the number of droplets. One can also add that this is the 
true value of k because namely this value corresponds to the sense of minimal 
work when we have the low dispersion. It corresponds to a minimal (in a 
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certain sense) entropy and later we can get additional work from increase of 
entropy. 

Dispersion at k — 0.6 will be 



Now we shall calculate the value of dispersion more accurate. Due to the 
similarity of nucleation the first cycle doesn't differ from the whole period. 
Function f3 for the first cycle will be 



Calculations for k x = 0.6 * 0.6 = 0.36 and k = 0.6 give (3 = (3' = 0.0255. 

In the most rough approximation one has to add f3' = 0.0255 to the 
previous value of f3 — 0.305 which leads to (3 = 0.32 -h 0.33. The smallness 
of /3' in comparison with f3 allows to use this linear approximation. For this 
value of (3 the value of dispersion will be 



It is interesting that this result can be gotten by another (rather artificial) 
procedure: 

We suppose that f3\ and fa are the first two terms of some series. We don't 
know other terms, but it is reasonable to suppose that the series resembles 
geometric progression with denominator fo/Pi- This leads to dispersion 



As it follows from fig. 19 the value of extremum remains k = 0.6. 

The absence of the shift of extremum is important and is necessary for 
this approach to be a self consistent. 

Now we shall see whether some other approaches can lead to essential 
reconsideration of result for dispersion. 

A way to make results more precise is take into account the shift of 
dispersion directly in initial formulas. Having written for the dispersion 
correction in the first cycle 



D°° = 2iV°°0.66 



(3 = (3i - (3 2 




D°° = D^ = 2iV°°0.64 



D°° = 2iy°°0.64i07 



£>oc = 2AToo, 



'7 
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Figure 19: Dispersion. 



with parameter 7, we get for the final distribution 

~ [°° dN l( lN 2 ...dN P n exp(- (iV ;~,f l)2 ) 



i=l 



cxp = — J — 

2(iV( fe ) - PN ± ) 

where P is the number of elementary intervals until the argument k, N( k "> is 
the number of droplets calculated on the base of the theory with averaged 
characteristics, Ni the mean number of droplets formed during interval num- 
ber % with account of fluctuations from previous intervals. The values a\ 
are given by 



j=p+i 
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and M is the total number of intervals. 

Having fulfilled integration Jf^ dN 1 dN 2 ...dN P , we get for a limit value of 
dispersion 

a a 

The approximate similarity of spectrums leads to equation on 7, which 
can be easily solved 



Calculations lead to 

7 (jfe = 0.6) = 0.51 

This result is very strange. It radically differs from the previous one. 
Certainly we made an error. The reason of the error in previous approach is 
that the duration of the first cycle is limited by k. So, we have to limit the 
duration of a whole period. The limit of the whole nucleation is, evidently, 
~ 1. 

The limit of integration corresponds to the current moment of time. It 
can not be greater than 1. So, we have to take it equal to 1. Then we have 
to recalculate (3 as 

^initial = Pi ~ @2 

P 1 = %l\ t dr{i-r feM-i A ) 

Jk JO 

fa = 16 f 1 d£ f 1 d V f dr(r - 3 (r - t?) 3 exp(-£ 4 ) exp( - V 4 ) 

Jk Jk Jo 

To calculate the value of dispersion we can act in two manners. The first 
way is to calculate dispersion at the point of extremum of f3. Now we shall 
show a dependence Pinitiai on k at fig. 20. 

Calculations give Pinitiai = 0.18 and for the final dispersion 

= 2iV oo 0.65 

This value practically coincides with a previous approach. 

It is necessary to stress that one can not directly use extremal properties 
of 7 to get p which provides /3 extremum. In reality the final characteristic 
is 7, but the calculations show that 7 has no extremum. Certainly, this is 
the weak point of approach based on extremal properties. 
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Figure 20: (^initial as a function of k 

Another way is to use k which provides the most precise value for the 
total number of droplets in this model. Having compared the mean total 
number of droplets in TAC calculated in this model one can see that this 
number has minimum 0.927 which is greater than the precise value and the 
value 0.90 given in the first iteration. This extremum is attained at k = 0.55. 
Namely this value will be chosen as k and this leads to 

Df = 0.62 

which lies in frames of precision of numerical simulation. 

One can add that there is no need to use k corresponding to extremum of 
the number of droplets calculated up to 1 instead of oo. Then the minimum 
is 0.852 and it is attained at k = 0.45. Actually there is no reasons to take 
this value because here the value of minimum is strongly less than the precise 
value and it corresponds to the maximum of deviation from the precise value. 
So, there are no reasons to take this value. 
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One has to stress that the procedure adopted here is really necessary. 
The two-period model with a fixed boundary in time scale used in jl] can 
not be the base of the correct calculation of dispersion. The reasons are the 
following 

• The behavior of supersaturation here is the model behavior (MB) with 
parameters, which characterize the influence of the first part of spec- 
trum. In investigation of stochastic effects the fluctuations of these pa- 
rameters are the source of fluctuations of the total number of droplets. 
The fuctional form of MB is chosen to ensure the correct number of 
droplets in frames of TAC. The change of parameters in MB leads to 
the imaginary change of external conditions in TAC which corresponds 
to the MB with given parameters. 

• The change of external parameters in TAC has to lead to the change 
of the boundary between parts. But this can not be done in the model 
with a fixed boundary. 

In our approach the summation of geometric progression and all other 
approaches ensure the possibility to overcome the restrictions of the two parts 
approach with a fixed boundary. The summation of geometric progression 
describes the equivalence of all points which is prescribed by the property of 
similarity of spectrums. 

The relative smallness of numerical errors in pP is caused by the following 
reasons 

• The moderate possibility of linearization (the non-linearity isn't too 
big) in kinetics of decay 

• The weak dependence N tot on p in frames of TAC 

• The weak dependence of N tot on N(p) in frames of TAC 

• The effects of stabilization are very strong. Really, instead of addition 
0.5 the benefit of the second part to dispersion is only 0.64 — 0.5 = 0.14 
(or even less if we use p ~ 0.55 instead of p ~ 0.5) 

• The existence of the special buffer part of nucleation period which will 
be described later 
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4.5 Numerical results 



Numerical simulation of nucleation can be done by the following method. We 
split the nucleation interval into many parts (steps). At every step a droplet 
will be formed or not. The probability to appear must be rather low, then we 
ensure the smallness of probability to have two droplets at the same interval. 
This means that the interval is "elementary". 

The process of formation is simulated by a random generator in a range 
[0, 1]. If a generated number is smaller than a threshold parameter u, then 
there will be no formation of a droplet. If it is greater than a threshold, we 
shall form a droplet. As a result we have spectrum / of droplets sizes. Now 
it is a chain of and 1. The parameter u descends according to macroscopic 
law [3] 

u = u exp(-TG/$) 

from a theory with averaged characteristics (it is based only on a conservation 
law without any averaging and can be used). Here 

$ is the initial supersaturation, F c is a free energy of critical embryos forma- 
tion, v c is a number of molecules inside a critical embryo, G is the number 
of molecules in a new phase taken in units of a molecules number density 
in a saturated vapor. By renormalization one can take away all parameters 
except G. 

To simplify calculations radically one can use the following representation 
for G: 

G — z^Gq — 3z G\ -\- 3z(j2 — C3 
where z is a coordinate of a front of spectrum, and Gi are given by 



Gi= [ f(x)x l dx 
Jo 



We needn't to recalculate G iy but can only ascend the region of integration, 
having added to integrals the functions z l f(z)dx at every step. 

Our results are given below. The interval is split into 30000 parts. Pa- 
rameter uo have been varied from up to 1 which leads to a different number 
of droplets. It is clear that the limit values are not good: at there are no 
droplets in the system, at 1 our intervals are not elementary. At every u 
results were averaged over 1000 attempts. 
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Figure 21: Shifts of droplets numbers as a function of In jV(°°) 

Shifts of droplets numbers are drawn at fig. 21 as a function of In jV(°°) 
It is seen that an analytical result about negligible value of corrections is 
correct. 

Dispersion as a function of In jV(°°) is shown in fig. 22. 

It is seen that the analytical value of dispersion coincides with numerical 
simulation. The ends of the curve correspond to a zero number of droplets 
and to a giant number of droplets when the elementary intervals are not 
elementary and have to be thrown out. 

Stochastic effects in dynamic conditions j3] can be analyzed by the same 
method. We needn't to describe it here. Numerical results are drawn below. 
Fig. 23 shows the shift in the number of droplets. It is small. Dispersion is 
drawn in fig. 24 (i.e. the value of 7). It is greater than in the case of decay. 

The physical reasons for the smallness of the droplet number shift for 
decay and for dynamic conditions will be different. 

For decay the reason is the following. The system wait the first droplet 
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Figure 22: Relative dispersion 7 as a function of In iV(°°) 

as long as necessary. Actually the time for kinetics of this system is t(G) 
with no connection with real time (certainly, the rate of nucleation has such 
connection). This phenomena is the reason for a smallness. 

In dynamic conditions there is a time dependent parameter - the intensity 
of external source. So, there is no such a reason. 

But here in the theory with averaged characteristics there is a property of 
a weak dependence of the total number of droplets on microscopic corrections 
for a free energy The same is valid also for fluctuation deviations. So 
there will be a weak effect of stochastic nucleation. 

Because the reasons for smallness of effect in decay and dynamic condi- 
tions are different it is interesting to see whether they continue to act when 
the supersaturation is stabilized at some moment. Analytical results shows 
that the will be an overlapping of two reasons. 

Really, if stabilization takes place at the period where the main consumers 
of vapor are going to appear then the majority of droplets appear in the 



56 



14i 




Figure 23: Shift in a droplets number as a function of lniV(°°) for dynamic 
conditions 

situation when there is no influence on the system. Then the situation for 
these droplets resembles decay conditions (and may be even better because 
the external supersaturation jT2] is going to decrease). So the reason for the 
decay situation works here. 

If stabilization takes place at the second cycle, then the behavior of su- 
persaturation is governed by droplets formed in dynamic conditions and we 
have here the reason for smallness in dynamic conditions. In both situations 
the effect is small. Numerical results confirm this conclusion. 

4.6 Conclusions 

The main result of this publication is a correct definition of all main charac- 
teristics of stochastic nucleation. It is shown that the main role in stochastic 
effects belongs to all droplets, but not to the main consumers of vapor. Only 
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Figure 24: Relative dispersion 7 as a function of lniV^ 00 ) for dynamic con- 
ditions 

the property of the nucleation conditions similarity allows us to solve the 
problem of account of all influences during the nucleation period. 

When all disadvantages of pQ, ^U] are shown it is clear that these publi- 
cations can not be considered as a solid base for nucleation investigation. 

But why results obtained in pQ, JU] are so close numerically to real values? 
The reason is that on a level of averaged characteristics there is a universality 
of nucleation process. So, the errors of PQ, cannot lead to a qualitatively 
wrong results. 

One has to stress that all corrections obtained in this paper are also 
universal ones. Certainly, they are some coefficients in decompositions and 
the functional form of decomposition is prescribed now. 

There is also a second specific reason for the smallness of an error in 
numerical values presented in PQ. The reason is the following 

• The process of nucleation can be split in three sequential parts. 
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• The first part is the part where the main consumers of vapor were 
appeared. Here the vapor depletion is small 

• The second (buffer) part is the part where conditions are ideal, the 
vapor depletion is small, but the droplets formed in this part can not 
attain big sizes even at the end of the whole nucleation period and, thus, 
can not consume enough vapor and they are not the main consumers 
of vapor 

• The third part is the part where the depletion isn't small and droplets 
appeared in this part aren't the main consumers of vapor in the nucle- 
ation period. 

• All parts have the lengths of one and the same scale. 

Certainly this structure was not declared in |T] which made the derivation in 
[T] illegal. 

The existence of the buffer part is necessary to balance the errors appeared 
from the fact that the fluctuations leads to to the absence of applicability of 
functional approximations for the nucleation kinetics. More correctly is to 
use the property of the "internal time of decay" which will be done separately. 

It seems that all effects considered here are negligible. For simple systems 
it is really true. But for systems with more complex kinetic behavior these 
effects can be giant. One of such systems is already described theoretically 
and this description will be presented soon in a separate publication. 

To compare results given here with the previous approach one can simply 
recall that the real error is the error in the relative deviation of dispersion 
from the standard value. In £Q this error is more than one quarter. Here the 
error is practically absent. 

In diffusion regime of droplets growth one has to use another approach 
based on jl], [T3]. In [T3] an explicit description of nucleation with account 
of stochastic effects was constructed. 
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